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Abstract 
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1 Introduction. 



It is well known that CP - violation in K° — K° mixing is described by the parameter 
e. Within the SM, this parameter is given by box diagrams. It depends in particular on 
the CKM matrix elements. On the other hand, the experimentally measured parameters 
are e and e'. e and e' enter the measured ratios of decay amplitudes of kaons into mi 
states. These amplitudes are superpositions of amplitudes A(K° — > (7111)1) = Aje l&1 of 
kaon decays into states with definite isospin I = 0, 2, Aj are weak amplitudes, 5i are 
strong rescattering phases(see Appendix for details). The parameter e can be expressed 
as pQ: 

e = e + % - . (1) 

Within the SM, ImAo originates from the so-called strong penguin diagrams. Ampli- 
tude A 2 also has an imaginary part which originates from electro-weak penguin diagrams. 

That is why ImAo >> ImAo. The ratio is much smaller than e and when the fit 

ReA 

of the CKM matrix parameters is performed, one equates the experimentally measured 
value of |e| and theoretical expression for |e|, neglecting the term — see |2],!3]- I n 

particular it was claimed in 4] that the contribution of is "at most a 2% correction 

l-! ReA 
to e". The aim of the present paper is to take this usually neglected term into account. 

In order to estimate the ratio — ^° we exploit the fact that it enters the expression 

e' ° 
for -IT]: 

e 



1 _ 1 c i(&2-5 ) 1 

e V2 e 



ImA 2 ImA 

w 

ReA ReA 



(2) 



ReA 2 

where w 



ReA 
e' 

The ratio — is experimentally measured and great amount of work was done in order to 
e J A 

calculate itfsee 13] - [TU] and refs. therein). In particular the quantity was computed 

lit . In 

theoretically using different methods. We shall use the results of these computations. 

We shall imply the following three step procedure for estimating — . 

ReA 

At first step we neglect ImA . Then \e t heor.\ coincides with \e exp , |, and we reproduce 
the results of 

At second step we take into account that ImAo 7^ 0, but neglect the contribution 

of EW penguins in Eq.(l2"l). Then we extract the value of from experimentally 

ReA 

t > 

measured quantity — with the help of Eq.(J2J). 

e 

At third step we take into account the contribution of EW penguins: ImA 2 7^ 0. The 

consequence is that one cannot extract from Eq.O). So one has to use the results 

ReA 
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of theoretical computation of 



ImAo 
ReA ' 



Finally, we perform a fit of CKM matrix parameters, taking the term 
into account and using numerical estimate of it, obtained at step 3. 



ImAo 
t~ in 

ReAn 



Eq.jl} 



2 Difference between e and e. 

The quantities e and e are related by Eq.((TJ). Taking into account that the phase of e is 
approximately | pQ (see also Appendix), from Eq.(JTJ we deduce: 



e + i 



ImAo 



ReA 



\2 



|e| 2 + 



1 ImA 
y/2 + ReA 



Thus: 



(3) 



1 ImAc 



+ 



\ 



, 2 _ 1 { ImAg 
C| 2 V.ReAo 



The experimentally measured value is [TTj: 



1 ImAo 
y/2 ReA ' 



(4) 



\e exp \ = 2.282(17) x 10" 3 . 



(5) 



Now we start our procedure of estimating |e|. At first step we neglect J 71 , and 



obtain: 



ReA 



c exp I 



2.282(17) x 10~ 3 . (6) 

This formula is always used in the fits of CKM matrix parameters, see 

Second step: We take into account that ImA ^ but neglect ImA 2 . Then Eq.(j2} 



reduces to: 



€ 
€ 



i - fx A 7r \ %u ImAo 
lei ReAn ' 



V2 



(7) 



for 



Taking into account that (So — $2)exp = 42 ±4° we obtain the following expression 
ImAo 

ReA ' 



ImAo 



V2\e\ e> 



ReA w e 

Substituting experimental values from [TT] we get: 



- = 1.8(4) x 1(T 3 , w = 0.045, lei = 2.282(17) x 10" 
e 

ImAo , \ a 

— —7 - = - 1.3 ±0.3 x 10~ 4 . 
ReAo 



(9) 



In this way we get the following value of |e|, which is the result of the second step: 
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|e| = 2.37(2) x 10~ 3 . (10) 

This number coincides with the value obtained in [13J, Eqs.(9.3),(9.4). 

Third step: Now let us take into account the presence of EW penguins: ImAi ^ 0. 

Then Eq.© does not allow to extract from the experimental data and we need 

ReA 

explicit theoretical result for — — — . As announced in the Introduction, such result was 

ReA 

obtained in the literature while calculating theoretically — . 

e 

e' 

In order to calculate — from Eq.(J2J) , one needs theoretical expressions for ImA^ and 

I1T1A2 (the values of ReA , ReA 2 , \e\, 5q — 62 and w are well measured experimentally). 

e' 

Short review of the history of — calculation can be found in [TO]- The expressions for 
ImAo and I171A2 are usually presented in the following form: 



ImA = -9Zl m (V td V t :)P^(l - Q IB) 



ImA 2 = -^Im{V u V: s )P^\ (11) 

where 

p {r) = Y,Vi(Qi)ii / = ' 2 - ( 12 ) 

% 

Here V t d and V t * are CKM matrix elements, Gp - Fermi constant, (Qi)o,2 are matrix 
elements of 4-quark operators responsible for K — > tttc decays, being their Wilson coef- 
ficients, ib introduces a correction due to isospin breaking effects: Qjb = —^—7 — t~—- 

w ImA 

From (fTT]) we have: 

^-T^™^' 01 ' 1 -^- (13) 

This formula contains the CKM matrix elements (which we are going to fit), but for 
the estimate of the small correction to |e| we can use mean values from J3]: Im(V t dV t * s ) = 
0.000127. ReA is well measured experimentally: ReA = 3.33 x 10~ 7 GeV. Concerning 

p(°\l — £Ijb), we use data from the calculations of — done in [S], which succeed in 



describing the experimental value of — . 

e 

Hadronic matrix elements were evaluated in 8\ using large N c - expansion. From Table 
2 of |H] we find the following range of values (corresponding to the quark condensate value 



(< >)' = 0.240 -0.260 GeV at /1 = 2 GeV): P^{1-Q IB ) = (7.1 ±2.1) x 10- 2 Ge1/ 3 . 

Substituting this into we get: = (-2.23 ± 0.66) x 10" 4 . 

RcAq 

This leads to the following range of values for |e|: 
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2.39 x 1(T 3 < e < 2.48 x 10~ 3 . (14) 

We have taken the paper [Hj as an example, and similar estimates can be made using 

e' 

other results, obtained in the framework of — calculation (see j5]-|10j). 
The range of values for |e| presented in Eq. fll4|) can be written as: 1 

|e| = (2.44 ±0.04) x 10" 3 , (15) 

and we use it in Section |3] to perform the fit of the parameters of CKM matrix. As we 
see the value of |e| is larger than that obtained at step 1 by (5 — 10)%. 



3 Fit of the parameters of CKM matrix 

We use in our fit of the CKM matrix experimentally measured values of modulus of matrix 
elements V ud ,V us ,V ub ,V cd ,V cs , V cb and also e, Am Bd and sin2(3. 

We assume these experimentally measured data to be normally distributed. Also the 
theoretical uncertainties are treated as normally distributed. Let us note that other people 
treat theoretical uncertainties in other way 

The most precise determination of \V u d\ comes from the averaging data from nuclear 
and neutron (3 decays [15] : 

\V ud \ = 0.9734 ±0.0008. (16) 

From kaon semileptonic decays the element |V^ S | is determined with the better accu- 
racy than in other methods (like hyperon semileptonic decays). We use the recent value 
H5|: 

\V US \ = 0.2196 ±0.0026. (17) 

From the inclusive and exclusive -B-decays governed by the transition b — > ul~vi we 
get [E]: 

\V uh \ = 0.0036 ±0.0007. (18) 

The element \V c d\ was measured in deep inelastic scattering of neutrinos and anti- 
neutrinos on nucleons with charm production [T5] : 

|V cd | = 0.224 ±0.016. (19) 

The best accuracy in |V^,| comes from the measurement of the ratio of hadronic W- 
decays to leptonic W-decays [15] : 

\V CS \ = 0.996 ±0.013. (20) 

The averaged value of \V c b\ extracted from exclusive and inclusive semileptonic B- 
decays including c quark is [T5] : 

IV^I = 0.041 ± 0.002. (21) 
1 We note that a number, very close to our central value, can be extracted from |14j . 
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Theoretical expression for |e| valid for m t > was first obtained in ^H]- I n modern 
notations it looks like: 

\i the °\ = ^^^ B ^ S ^ Xc)Im[(V cs V; d ) 2 ] + rj tt S(x t , x t )Im[{V ts V t * d f\ 
+27 lct S(x c , x t )Im[V cs V; d V ts V t * d }). (22) 
Here, the S(xi,Xj) are usually called the Inami-Lim functions [T?j : 
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(23) 



4(l-x i )(l- 

where Xj = mf/m^ depend on the masses of c quark and t quark (m c = 1.2±0.2 GeV 
mt = 174.3 ± 5.1 GeV [EJ, = 80.42 ± 0.04 GeV HE]). The QCD corrections have 
been calculated to next-to-leading order: 7] cc = 1.32 ± 0.32 [TH] , ?7» = 0.574 ± 0.01 [Uj, 
?7 rf = 0.47 ± 0.04 [21]. The kaon decay constant extracted from the K + — > |U + z/ decay 
width equals: fx = 160.4 ± 1.9 MeV [TH]. The K$ — Kl mass difference is Arrix = 
(3.491 ± 0.006) x 10~ 15 GeV |18j . The world average for the bag parameter B^ reads: 
B K = 0.87 ± 0.06 ± O.U quench [22\. Fermi constant G F = 1.16639(1) x W~ 5 GeV- 2 \T^. 

From the study of B° — B° oscillations the experimental value of \Vu\ should be 
extracted: 

Am B d = ^VBm Bd m 2 w S{x t )flB Bd \V td V^\\ (24) 

where rj B = 0.55 ± 0.01 [20] is a QCD correction, m Bd = 5.2794 ± 0.0005 GeV [H] is 
the B meson mass, mw is the W boson mass, S(xt) is the Inami-Lim function for the 

2 

box diagram, x t = — »jf-, fs d is the B meson decay constant, and Bs d is the so-called bag 

factor. We use the following numerical value: f Bd \fB~ d = 230 ± 28 ± 28 MeV [23]. 

From B decays to CP eigenstates containing charmonium and neutral K-meson sin 2(3 
is measured with good accuracy. The average result of Belle and BaBar is [24J: 

sin 2(3 = 0.73 ± 0.05(stat) ± 0.035(syst). (25) 

Theoretical formula for sin 2(3 comes from the consideration of the unitarity triangle: 

sm ^ = ^+(i-- P r (26) 

The x 2 expression which we minimize looks like: 



X 2 (A,\p, V )= p!* 22- + — 22- + 




Tjtheo Trexp \ 2 / irtheo \f ex P\ 2 

* us ' us 1 i / ufe ^ ub I 



cs cs I i / cb cb 




2 i\rtheo \/exp\ 2 



°Vc6 / 

' Amg l d eo - Am^f \ 2 ^ (e^ 60 ) - \e exp \ \ 2 , ( sin2(3 theo - sin2(3 exp ^ 2 



(27) 
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where theoretical expressions depend on the Wolfenstein parameters A, A, p, r\. Expression 
(|27|) was minimized varying A, A, p, 77. 

Performing the fit we use the value of |e* heo | from Eq.(|15|). The main uncertainty in 
originates from that in Bk and it dominates in a- t . That is why we use = 0.4 x 10~ 3 . 
Here are our results: 

A = 0.2229 ±0.0021 
A = 0.83 ±0.04 



~theo 



V 

p 



0.35! 



[-0.05 
-0.04 



20 +0 - 08 

u - zu -0.09 



X 2 /n.d.o.f. = 8.1/5 . 

In Fig.l you can see a set of bounds on the parameters p and fj of the CKM matrix. 
They comprise three circles, two branches of a hyperbola, and two straight lines. Three 
circles originate from the V u b measurement (the green one), the measurement of Am^ (the 
red one) and from the lower bound on Am# s (the yellow one). The hyperbola originates 
from the measurement of CP violation in the mixing of i^-mesons. Straight lines come 
from the measurement of CP asymmetry in B^(B^) — > J/^/K decays. 




Figure 1: The domains at (p,i]) plane allowed at lcr from V u b, Am^, Ek and sin 2(3 
measurements. 9ffl$CL. upper bound from the search of Am^ is shown as well. 



4 Conclusions 

Numerical difference of the quantities e and e (which describe CP- violation in K- mesons) 
was estimated in the Standard Model. Fit of CKM matrix patameters accounted for this 
difference was performed. 
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A Basic formulas for K°-K° system 

It is known that states K° and K° are not mass eigenstates. Mass eigenstates are their 
linear combinations: 



K 4 



1 



1+ I e I 2 
1 



V2 



V2 



1+ e 



- 12 



K° -K° K° + K° 



V2 



V2 



(28) 



Let's denote matrix elements of the effective Hamiltonian between K° and K° states 
as follows: 



< \ H \ >=< K \ H \ K >= M 



< K° I H I K° >= M 



12 



12) 



< K° | H | K° >= M* 2 - l -T\ 2 . 
The eigenvalues and eigenvectors of this matrix Hamiltonian are: 



(29) 



A± = M - l -v ± J (M 12 - ^r 12 )(Mf 2 - l -rt 2 ), 



IV1 12 2 i 12 



f M+ = pM° + qM° q 

\ M_ = pM° -qM° ' p~\M 12 - fr 12 ' 

Introducing quantity e according to the following definition: 

q 1 — e 
V 



we come to Eq. 

Taking into account that T 12 is real and 
expression: 



l + e 

ImM 12 



(30) 



ReM 12 



(31) 



0.1 [13^ we get the following 



q UmM\ 2 
-ml 

V 



Mi 



12 



|Ti 2 



(32) 



S 



Eigenvalues of Hamiltonian may be written as A± = (m± — |T ± ) 2 , where m± are 
masses of corresponding states and - their widths. Then denoting K + and K_ states 
as Ks and Kl respectively, we have A_ — A + = 2mK(mL — m s ~ 2 (^ L ~ ^s))- On the 



other hand A_ - A + = -2yJ(M 12 - fr 12 )(M 1 * 2 - §r? 2 ) ~ -2(M 12 - %T 12 ). This leads to: 

HmMi2/m K 



l-2e. 



p ( mL _ m5 _i(r L _r 5 )) 

Taking into account that Ts « Tl and Aijils ~ we obtain: 



(33) 



HmMi2/2mK 



~ a 1 4 ' 



(34) 



(Am L5 -f(r L -r 5 )) 

Thus calculating ImMi 2 within the SM we find the theoretical prediction for e. (Let 
us note that since ImM 12 is negative, the phase of e approximately equals |). 

Now we proceed to decays of kaons into pairs of pions, whose amplitudes are well 
measured experimentally. 

It is convenient to deal with the amplitudes of the decays into the states with definite 
isospin: 



A(K° - 


-> 7T + 7T ) 


A(K° - 


-> 7T + 7T~) 


A(K° - 


-> 71 71 ) 


A(K° - 


-> 71 7T J 
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(35) 



where "2" and "0" are the values of (tttt) isospin, £ 2 ,o are the (small) weak phases 
which originate from CKM matrix and 5 2 are the strong phases of 7T7T-rescattering. 
Experimentally measured quantities are: 



V+- = 
Voo = 



A{K L -> 7T + 7T-) 

A{K L - 



* 73 r 71 
7T 7T ) 



7r 7r ' 



A(X 5 

For the amplitudes of and Kg decays into tx+tx^ we obtain: 



(36) 



A{Kl _» 7r+7r - ) = 



^e i52 2i sin £ 2 + -^L v ^e Wb 2i sin £ 
v 3 v 3 



+ 



+ 



v/2 



^=e l<52 2 cos 6 + ^=V2e i5o 2 cos^o 
v 3 v 3 



A{Ks _» 7r+7r - ) = 



^=e i52 2cos6 + ^W2e Wo 2cos£ 
V3 v 7 ^ 



(37) 



where in the last equation we omit the terms which are proportional to the product of 
two small factors, e and sin£ ,2- For the ratio of these amplitudes we get: 



V+- 



A{K L 



IX 1 7T 



A(K S 



.sin£ i^ &2 ^ a 2 cos£ 2 
e + i — + 



7T 7T 



COs£, 







y/2 a cos £ 



sin £ 2 sin £ 



cos £ 2 cos £ 



where we neglect the terms of the order of (a 2 /a ) 2 sin£ 0i2 , because a 2 /a ~ 1/22. 
The analogous treatment of K L S — ► vr°7r° decay amplitudes leads to: 



Voo 



A (K L -> A ) = l + . sin^o _ I . c i( fa -j )^ fl2cos6 
A(Ks — ► 7r°vr°) cos£o aocos£o 



sin ^ 2 _ sin^o 
cos £ 2 cos £ 



Introducing conventional quantities e 



+ ±7700 and e' = \r] + ^ - \r] 00 , we get: 



1 



ItuAq 



[ImA 2 — w/mA ] 



where A 2 = e^ 2,0 a 2 0, and u> 



(38) 



«2 



Equations (jBEJ) are our starting point in the present paper; see Introduction. 
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